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Abstract 

We review the main tools which allow for the statistical characterization of 
weighted networks. We then present two case studies, the airline connection network 
and the scientific collaboration network, which are representative of critical infras- 
tructures and social systems, respectively. The main empirical results are (i) the 
broad distributions of various quantities and (ii) the existence of weight-topology 
correlations. These measurements show that weights are relevant and that in general 
the modeling of complex networks must go beyond topology. We review a model 
which provides an explanation for the features observed in several real-world net- 
works. This model of weighted network formation relies on the dynamical coupling 
between topology and weights, considering the rearrangement of weights when new 
links are introduced in the system. 

PACS numbers: 89.75.-k, -87.23.Ge, 05.40.-a 



1 Introduction 

Networked structures arise in a wide array of different contexts such as tech- 
nological and transportation infrastructures, social phenomena, and biological 
systems. These highly interconnected systems have recently been the focus of 
a great deal of attention that has uncovered and characterized their topo- 
logical complexity [1,2,3,4]. Along with a complex topological structure, real 
networks display a large heterogeneity in the capacity and intensity of the 
connections — the weight of the links. For example, in ecology the diversity of 
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the predator-prey interaction is believed to be a critical ingredient of ecosys- 
tems stability [5,6], and in social systems, the weight of interactions is very 
important in the characterization of the corresponding networks [7] . Similarly, 
the Internet traffic [3] or the number of passengers in the airline network [4,8,9] 
are crucial quantities in the study of these systems. 

In this paper we review a set of metrics combining weighted and topolog- 
ical observables that allows to characterize the complex statistical proper- 
ties of the strength of edges and vertices and to investigate the correlations 
among weighted quantities and the underlying topological structure. Specifi- 
cally, we present results on the scientific collaboration network and the world- 
wide air-transportation network, which are representative examples of social 
and large infrastructure systems, respectively. The measures on weighted net- 
works [9,10,11,12] have shown that they can exhibit additional complex prop- 
erties such as broad distributions and non-trivial correlations of weights that 
do not find an explanation just in terms of the underlying topological struc- 
ture. The heterogeneity in the intensity of connections may thus be very im- 
portant in real-world systems and cannot be overlooked in their description. 
Motivated by these observations, we review also a model for weighted networks 
we have recently proposed [13], which naturally produces topology- weight cor- 
relations and broad distributions for various quantities. 

2 Tools for the characterization of weighted networks 

We briefly review the different tools which allow for a first statistical charac- 
terization of weighted complex networks. 

• Weights 

The properties of a graph can be expressed via its adjacency matrix a^-, 
whose elements take the value 1 if an edge connects the vertex % to the 
vertex j, and otherwise (with i,j = 1,...,N, where N is the size of the 
network). Weighted networks are usually described by a matrix Wij speci- 
fying the weight on the edge connecting the vertices i and j (wij = if the 
nodes i and j are not connected). In the following we will consider only the 
case of symmetric positive weights ^ = Wji > 0. 

• Connectivity and weight distributions 

The standard topological characterization of networks is obtained by 
the analysis of the probability distribution P{k) that a vertex has de- 
gree k. Complex networks often exhibits a power-law degree distribution 
P{k) ~ A;" 7 with 2 < 7 < 3. Similarly, a first characterization of weights is 
obtained by the distribution P{w) that any given edge has weight w. 

• Weighted connectivity: Strength 

Along with the degree of a node, a very significative measure of the net- 
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work properties in terms of the actual weights is obtained by looking at the 
vertex strength Sj defined as [14,9,15] 

s i= (!) 
jev(i) 

where the sum runs over the set V(i) of neighbors of i. The strength of 
a node integrates the information both about its connectivity and the im- 
portance of the weights of its links, and can be considered as the natural 
generalization of the connectivity. When the weights are independent from 
the topology, we obtain that the strength of the vertices of degree k is 
s(k) ~ (w)k where (w) is the average weight. In the presence of correlations 
we obtain in general s(k) ~ Ak 13 with (3 = 1 and A ^ (w) or (3 > 1. 



• Weighted clustering 

The clustering coefficient q measures the local cohesiveness and is de- 
fined for any vertex % as the fraction of connected neighbors of % [16]. The 
average clustering coefficient C = N' 1 J2i Q thus expresses the statistical 
level of cohesiveness measuring the global density of interconnected vertex 
triplets in the network. Further information can be gathered by inspecting 
the average clustering coefficient C(k) restricted to the class of vertices with 
degree k. The topological clustering, however, does not take into account the 
fact that some neighbors are more important than others. In order to solve 
this incongruity we introduce a measure of the clustering that combines the 
topological information with the weight distribution of the network. The 
weighted clustering coefficient is defined as [9] 

W = ~7l TT 2- — o a i3 a ih^jh- (2) 

This quantity c w (i) counts for each triple formed in the neighborhood of 
the vertex % the weight of the two participating edges of the vertex %. In 
this way we are not just considering the number of closed triangles in the 
neighborhood of a vertex but also their total relative weight with respect 
to the vertex' strength. The factor s,i(k,i — 1) is a normalization factor and 
ensures that < cf < 1. Consistently, the cf definition recovers the topo- 
logical clustering coefficient in the case that Wij = const. It is customary to 
define C w and C w (k) as the weighted clustering coefficient averaged over all 
vertices of the network and over all vertices with degree k, respectively. In 
the case of a large randomized network (lack of correlations) it is easy to see 
that C w = C and C w (k) = C(k). In real weighted networks, however, we 
can face two opposite cases. If C w > C, we are in presence of a network in 
which the interconnected triples are more likely formed by the edges with 
larger weights. On the contrary, C w < C signals a network in which the 
topological clustering is generated by edges with low weight. In this case it 
is obvious that the clustering has a minor effect in the organization of the 
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network since the largest part of the interactions (traffic, frequency of the 
relations, etc.) is occurring on edges not belonging to interconnected triples. 
The same may happen for C w (k), for which it is also possible to analyze 
the variations with respect to the degree class k. 

• Weighted assortativity: Affinity 

Another quantity used to probe the networks' architecture is the average 
degree of nearest neighbors of a vertex i 

k n n,i = ~T ^3 ' (^) 

* jev(») 

where the sum runs on the nearest neighbors vertices of each vertex %. From 
this quantity a convenient measure to investigate the behavior of the degree 
correlation function is obtained by the average degree of the nearest neigh- 
bors, k nn (k), for vertices of degree A; [17]. In the presence of correlations, the 
behavior of k nn (k) identifies two general classes of networks. If k nn (k) is an 
increasing function of k, vertices with high degree have a larger probabil- 
ity to be connected with large degree vertices. This property is referred in 
physics and social sciences as assortative mixing [18]. On the contrary, a de- 
creasing behavior of k nn (k) defines dis assortative mixing, in the sense that 
high degree vertices have a majority of neighbors with low degree, while the 
opposite holds for low degree vertices. 

In the case of weighted networks an appropriate characterization of the 
assortative behavior is obtained by the weighted average nearest neighbors 
degree, defined as 

1 N 

k nn,i = —Y, a i3 W H k 3- ( 4 ) 
bt 3=1 

In this case, we perform a local weighted average of the nearest neighbor 
degree according to the normalized weight of the connecting edges, w^/sj. 
This definition implies that k™ n i > k nn ^ if the edges with the larger weights 
are pointing to the neighbors with larger degree and k™ n i < k nn ,i in the 
opposite case. The k™ n i function thus measures the effective affinity to con- 
nect with high or low degree neighbors according to the magnitude of the 
actual interactions. As well, the behavior of the function k™ n (k) (defined 
as the average of k™ ni over all vertices with degree k), marks the weighted 
assortative or disassortative properties considering the actual interactions 
among the system's elements. 

• Disparity 

For a given node % with connectivity fcj and strength Sj different situations 
can arise. All weights Wij can be of the same order Sj/fcj. In contrast, the 
most heterogeneous situation is obtained when one weight dominates over 
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all the others. A simple way to measure this "disparity" is given by the 
quantity Y 2 introduced in other context [19,20] 



^ 2 «= E 

iev(i) 



(5) 



If all weights are of the same order then Y 2 ~ (for fc, ^> 1) and if a 
small number of weights dominate then Y 2 is of the order 1/n with n of 
order unity. This quantity was recently used for metabolic networks [21] 
which showed that for these networks one can identify dominant reactions. 



3 Empirical results 

3.1 Weighted networks data 

Prototypical examples of weighted networks can be found in the world-wide 
airport network (WAN) [8,9] and the scientific collaboration network (SCN) [22,23]. 
In the airport network each given weight Wij is the number of available seats 
on direct flight connections between the airports % and j. For the WAN, we 
analyze the International Air Transportation Association (IATA) 1 database 
containing the world list of airports pairs connected by direct flights and the 
number of available seats on any given connection for the year 2002. The result- 
ing air-transportation graph comprises N = 3880 vertices denoting airports 
and E = 18810 edges accounting for the presence of a direct flight connection. 
The average degree of the network is (k) = 2E/N = 9.70, while the maximal 
degree is 318. 

In the SCN the nodes are identified with authors and the weight depends on 
the number of co-authored papers [22,9]. We consider the network of scien- 
tists who have authored manuscripts submitted to the e-print archive relative 
to condensed matter physics (|http: / /xxx.lanl.gov/ archive/cond-mat ) between 



1995 and 1998. Scientists are identified with nodes and an edge exists be- 
tween two scientists if they have co-authored at least one paper. The resulting 
connected network has iV = 12722 nodes, with an average degree (i.e. av- 
erage number of collaborators) (k) = 6.28 and maximal degree 97. For the 
SCN we follow the definition of weight introduced in Ref. [22]: The inten- 
sity Wy of the interaction between two collaborators i and j is defined as 
= J2 P ^i^/( n p ~ 1) where the index p runs over all papers, n p is the 
number of authors of the paper p, and 5f is 1 if author i has contributed to 
paper p, and otherwise. This definition seems to be rather objective and 
representative of the scientific interaction: It is large for collaborators having 
many papers in common but the contribution to the weight introduced by any 
given p aper is inversely proportional to the number of authors. 
1 http://www.iata.org. 
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3.2 Empirical results 

3. 2. 1 Topological properties 

The topological properties of the SCN network and other similar networks 
of scientific collaborations have been studied in Ref. [22] and we report on 
Fig. (1A) the connectivity distribution showing a relatively broad law. 

As shown in Fig. (IB), the topology of the WAN exhibits both small-world and 
scale-free properties as already observed in different dataset analyses [10,8]. In 
particular, the average shortest path length, measured as the average number 
of edges separating any two nodes in the network, shows the value (£) = 4.37, 
very small compared to the network size N. The degree distribution, on the 
other hand, takes the form P(k) = k~ n f(k/k x ), where 7 ~ 2.0 and f(k/k x ) is 
an exponential cut-off function that finds its origin in physical constraints on 
the maximum number of connections that a single airport can handle [4,8]. 
The airport connection graph is therefore a clear example of heterogeneous 
network showing scale-free properties on a definite range of degree values. 

3.2.2 Strength distribution 

The probability distribution P(s) that a vertex has strength s is heavy tailed 
in both networks and the functional behavior exhibits similarities with the 
degree distribution P(k) (see Fig. 1). A precise functional description of the 
heavy-tailed distributions may be very important in understanding the net- 
work evolution and will be deferred to future analysis. This behavior is not 
unexpected since it is plausible that the strength Sj increases with the vertex 
degree ki, and thus the slow decaying tail of P(s) stems directly from the very 
slow decay of the degree distribution. 

3.2.3 Topology-weight correlations 

In Fig. 2 we report the behavior obtained for both the real weighted net- 
works and their randomized versions, generated by a random re-distribution 
of the actual weights on the existing topology of the network. For the SCN 
the curves are very similar and well fitted by the uncorrelated approximation 
s(k) = (w)k. Strikingly, this is not the case of the WAN. Fig. 2B clearly shows 
a very different behavior for the real data set and its randomized version. In 
particular, the power-law fit for the real data gives an "anomalous" exponent 
/Swan = 1.5 ± 0.1. This implies that the strength of vertices grows faster than 
their degree, i.e. the weight of edges belonging to highly connected vertices 
tends to have a value higher than the one corresponding to a random assign- 
ment of weights. This denotes a strong correlation between the weight and the 
topological properties in the WAN, where the larger is an airport, the more 
traffic it can handle. 
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Fig. 1. A Degree and strength distribution in the scientific collaboration net- 
work. The degree k corresponds to the number of co-authors of each scientist 
and the strength represent its total number of publications. The distributions are 
heavy-tailed even if it is not possible to distinguish a definite functional form. B 
The same distributions for the world-wide airport network. The degree is the num- 
ber of non-stop connections to other airports and the strength is the total number 
of passengers handled by any given airport. In this case, the degree distribution 
can be approximated by the power-law behavior P{k) ~ fc -7 with 7 = 1.8 =L 0.2. 
The strength distribution has a heavy-tail extending over more than four orders of 
magnitude. 



The fingerprint of these correlations is also observed in the dependence of the 
weight Wij on the degrees of the end point nodes fc* and kj. For the WAN 
the behavior of the average weight as a function of the end points degrees 
can be well approximated by a power-law dependence ~ {kikj) 6 with 

an exponent 9 = 0.5 ± 0.1. This exponent can be related to the j3 exponent 
by noticing that s(k) ~ k(kkj) 6 , resulting in (5 = 1 + 9, if the topological 
correlations between the degree of connected vertices can be neglected. This is 
indeed the case of the WAN where the above scaling relation is well satisfied 
by the numerical values provided by the independent measurements of the 
exponents. In the SCN, instead, is almost constant for over two decades 
confirming a general lack of correlations between the weights and the vertices 
degree. In this case 9 = and the relation (3 = 1 + 9 also holds. 
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Fig. 2. Average strength s(k) as function of the degree k of nodes. A In the 
scientific collaboration network the real data are very similar to those obtained in a 
randomized weighted network. Only at very large k values it is possible to observe a 
slight departure from the expected linear behavior. B In the world airport network 
real data follow a power-law behavior with exponent (3 = 1.5 ± 0.1. This denotes 
anomalous correlations between the traffic handled by an airport and the number 
of its connections. 

3.2.4 Weighted clustering and assortativity 

We present the results [9] obtained for both the SCN (see Fig.3) and the WAN 
(see Fig. 4) by comparing the regular topological quantities with the weighted 
ones introduced above. In the figures we report the relative difference of the 
values obtained for the topological and weighted quantities. It is striking to 
observe that for large degree values we observe up to 100% relative difference 
signaling a strong difference in the clustering and correlation properties if we 
take into account the weighted nature of the networks. 

• SCN 

■ (i) The measurements indicate that the SCN has a monotonously decay- 
ing spectrum C(k). This implies that hubs present a much lower clus- 
tered neighborhood than low degree vertices which can be interpreted as 
the evidence that authors with few collaborators usually work within a 
well defined research group in which all the scientists collaborate together 
(high clustering). Authors with a large degree, however, collaborate with 
different groups and communities which on their turn do not have often 
collaborations, thus creating a lower clustering coefficient. 
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Fig. 3. Comparison of topological and weighted quantities for the SCN. A The 
weighted clustering separates from the topological one around k > 10. This marks 
a difference for authors with larger number of collaborators. B The assortative 
behavior is enhanced in the weighted definition of the average nearest neighbors 
degree. It is worth remarking the large relative difference (up to 50-100%) of the 
weighted and topological quantities (lower part of the figures). 



• (ii) The inspection of C w (k) shows generally that for k > 10 the weighted 
clustering coefficient is larger than the topological one. This implies that 
authors with many collaborators tend to publish more papers with inter- 
connected groups of co-authors and is a signature of the fact that influen- 
tial scientists form stable research groups where the largest part of their 
production is obtained. 

• (iii) Furthermore, the SCN exhibits an assortative behavior in agreement 
with the general evidence that social networks are usually denoted by a 
strong assortative character [18]. Finally, the assortative properties find 
a clearcut confirmation in the weighted analysis with a k™ n (k) strikingly 
growing as a power-law function of k. 

• WAN 

A different picture is found in the WAN, where the weighted analysis 
provides a richer and somehow different scenario. 

• (i) This network also shows a decaying C(k), consequence of the role of 
large airports that provide non-stop connections to very far destinations on 
an international and intercontinental scale. These destinations are usually 
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Fig. 4. Topological and weighted quantities for the WAN. A The weighted cluster- 
ing coefficient is larger than the topological one in the whole degree spectrum. B 
knn(k) is reaching a plateau for k > 10 denoting the absence of marked topological 
correlations. On the contrary k™ n (k) exhibits a more definite assortative behavior. 
Also in this case the relative difference of the weighted and topological quantities 
are of the order of 100% (lower part of the figures). 

not interconnected among them, giving rise to a low clustering coefficient 
for the hubs. 

• (ii) We find, however, that C w /C ~ 1.1, indicating an accumulation of 
traffic on interconnected groups of vertices. 

• (iii) The weighted clustering coefficient C w (k) has much more limited 
variation in the whole spectrum of k. This implies that high degree air- 
ports have a progressive tendency to form interconnected groups with high 
traffic links, thus balancing the reduced topological clustering. Since high 
traffic is associated to hubs, we have a network in which high degree nodes 
tend to form cliques with nodes with equal or higher degree, the so-called 
rich-club phenomenon [24]. 

• (iv) The topological k nn (k) does show an assortative behavior only at 
small degrees. For k > 10, k nn (k) approaches a constant value, a fact 
revealing an uncorrelated structure in which vertices with very different 
degrees have a very similar neighborhood. The analysis of the weighted 
k™ n (k), however, exhibits a pronounced assortative behavior in the whole 
k spectrum, providing a different picture in which high degree airports 
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have a larger affinity for other large airports where the major part of the 
traffic is directed. 



4 Modeling weighted networks 

Previous approaches to the modeling of weighted networks focused on growing 
topologies where weights were assigned statically, i.e. once and for ever, with 
different rules related to the underlying topology [14,25] (see also [26] for 
a more recent model). These mechanisms, however, overlook the dynamical 
evolution of weights according to the topological variations. We can illustrate 
this point in the case of the airline network. If a new airline connection is 
created between two airports it will generally provoke a modification of the 
existing traffic of both airports. In general, it will increase the traffic activity 
depending on the specific nature of the network and on the local dynamics. In 
the following, we review a model that takes into account the coupled evolution 
in time of topology and weights. Instead of drawing randomly the weights, 
an alternative consists in coupling the evolution of the weights and of the 
topology and allowing the dynamical evolution of weights during the growth 
of the system. This mimics the evolution and reinforcements of interactions 
in natural and infrastructure networks. 

The model dynamics starts from an initial seed of Nq vertices connected by 
links with assigned weight wo- At each time step, a new vertex n is added with 
m edges (with initial weight w ) that are randomly attached to a previously 
existing vertex i according to the probability distribution 

n_ = (6) 

This rule of "busy get busier" relaxes the usual degree preferential attachment, 
focusing on a strength driven attachment in which new vertices connect more 
likely to vertices handling larger weights and which are more central in terms of 
the strength of interactions. This weight driven attachment (Eq. (6)) appears 
to be a plausible mechanism in many networks (see also [27] for an alternative 
mechanism in which edges with large weight are chosen for preferential attach- 
ment). In the Internet new routers connect to more central routers in terms of 
bandwidth and traffic handling capabilities and in the airport networks new 
connections are generally established to airports with a large passenger traffic. 
Even in the SCN this mechanism might play a role since an author with more 
co-authored papers is more visible and open to further collaborations. 

The presence of the new edge (n, i) will introduce variations of the existing 
weights across the network. In particular, we consider the local rearrangements 
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Fig. 5. Illustration of the construction rule. A new node n connects to a node % 
with probability proportional to s, / J2j s j ■ The weight of the new edge is wq and the 
total weight on the existing edges connected to i is modified by an amount equal to 
S. 

of weights between % and its neighbors j G V(i) according to the simple rule 
w i:j -> Wij + Awy, (7) 



where 



Awij = 8^-. (8) 

Si 



This rule considers that the establishment of a new edge of weight wq with the 
vertex % induces a total increase of traffic S that is proportionally distributed 
among the edges departing from the vertex according to their weights (see 
Fig. 5), yielding Sj — > s$ + 5 + w . At this stage, it is worth remarking that 
while we will focus on the simplest model with 5 = const, different choices 
of Awij with heterogeneous 5i or depending on the specific properties of each 
vertex ki, Sj) can be considered [28,29,30]. Finally, after the weights have 
been updated the growth process is iterated by introducing a new vertex with 
the corresponding re-arrangement of weights. 

The model depends only on the dimensionless parameter 5 (rescaled by wo), 
that is the fraction of weight which is 'induced' by the new edge onto the oth- 
ers. According to the value of S, different scenarios are possible. If the induced 
weight is 5 ~ 1 we mimic situations in which an appreciable fraction of traffic 
generated by the new connection will be dispatched in the already existing 
connections. This is plausible in the airport networks where the transit traffic 
is rather relevant in hubs. In the case of 5 < 1 we face situations such as the 
SCN where it is reasonable to consider that the birth of a new collaboration 
(co-authorship) is not triggering a more intense activity on previous collabo- 
rations. Finally, 5 > 1 is an extreme case in which a new edge generates a sort 
of multiplicative effect that is bursting the weight or traffic on neighbors. 

The network's evolution can be inspected analytically by studying the time 
evolution of the average value of Si(t) and ki(t) of the i-th vertex at time 
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t, and by relying on the continuous approximation that treats k, s and the 
time t as continuous variables [1,2,13]. The behavior of the strength and the 
connectivity are easily obtained and one has in the long time limit 

Si (t)~(26 + l)ki(t) (9) 



This proportionality relation s ~ k implies (3 = 1 but the prefactor is differ- 
ent from (w) which indicates the existence of correlations between topology 
and weights. This relation (9) is particularly relevant since it states that the 
weight-driven dynamics generates in Eq. (6) an effective degree preferential at- 
tachment that is parameter independent. This highlights an alternative micro- 
scopic mechanism accounting for the presence of the preferential attachment 
dynamics in growing networks. 

The behavior of the various statistical distribution can be easily computed 
and one obtains in the large time limit P(k) ~ A; -7 and P(s) ~ s~ 7 with 

45 + 3 

1 = WTv < 10 > 



This result shows that the obtained graph is a scale-free network described 
by an exponent 7 G [2, 3] that depends on the value of the parameter 5. In 
particular, when the addition of a new edge doesn't affect the existing weights 
(5 = 0), the model is topologically equivalent to the Barabasi-Albert model 
[31] and the value 7 = 3 is recovered. For larger values of 5 the distribution 
is progressively broader with 7^2 when 5 — > 00. This indicates that the 
weight-driven growth generates scale-free networks with exponents varying in 
the range of values usually observed in the empirical analysis of networked 
structures [1,2,3]. Noticeably the exponents are non- universal and depend 
only on the parameter 5 governing the microscopic dynamics of weights. The 
model therefore proposes a general mechanism for the occurrence of varying 
power-law behaviors without resorting on more complicate topological rules 
and variations of the basic preferential attachment mechanism. 

Similarly to the previous quantities, it is possible to obtain analytical expres- 
sions for the evolution of weights and the relative statistical distribution [13]. 
The probability distribution P(w) is in this case also a power-law P(w) ~ w~ a 
where a = 2 + |. The exponent a has large variations as a function of the 
parameter 5 and P(w) moves from a delta function for S = to a very slow 
decaying power-law with a = 2 if 5 — > 00. This feature clearly shows that 
the weight distribution is extremely sensitive to changes in the microscopic 
dynamics ruling the network's growth. 

In summary, the networks generated by the model display power-law behavior 
for the weight, degree and strength distributions with non-trivial exponents 
depending on the unique parameter defining the model's dynamics. These 
results suggest that the inclusion of weights in networks modeling naturally 
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explains the diversity of scale-free behavior empirically observed in real net- 
worked structures. Strikingly, the weight-driven growth recovers an effective 
preferential attachment for the topological properties, providing a microscopic 
explanation for the ubiquitous presence of this mechanism. 

5 Conclusions and perspectives 

A more complete view of complex networks is thus provided by the study of 
the interactions defining the links of these systems. The weights characterizing 
the various connections exhibit complex statistical features with highly vary- 
ing distributions and power-law behavior. In particular we have considered 
the specific examples of the scientific collaboration and world-wide airport 
networks where it is possible to appreciate the importance of the correlations 
between weights and topology in the characterization of real networks prop- 
erties. Indeed, the analysis of the weighted quantities and the study of the 
correlations between weights and topology provide a complementary perspec- 
tive on the structural organization of the network that might be undetected 
by quantities based only on topological information. The weighted quantities 
thus offer a quantitative and general approach to understanding the com- 
plex architecture of real weighted networks. The empirical results shows that 
purely topological models are inadequate and that there is a need for mod- 
els going beyond pure topology. The model we have presented is possibly the 
simplest one in the class of weight-driven growing networks. A novel feature 
in the model is the weight dynamical evolution occurring when new vertices 
and edges are introduced in the system. This simple mechanism produces a 
wide variety of complex and scale-free behavior depending on the physical 
parameter 5 that controls the local microscopic dynamics. While a constant 
parameter 5 is enough to produce a wealth of interesting network properties, 
a natural generalization of the model consists in considering S as a function of 
the vertices degree or strength. Similarly, more complicated variations of the 
microscopic rules may be implemented to mimic in a detailed fashion particu- 
lar networked systems [28,29,30,32]. In particular, an important feature to be 
included in the description of the airline network is the geographical embed- 
ding of the network [33]. In this perspective the present model appears as a 
general starting point for the realistic modeling of complex weighted networks. 
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